Fermi edge polaritons in a highly degenerate 2D electron gas: a diagrammatic theory 
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We present a theoretical study on polaritons in highly doped semiconductor microcavities. In 
particular, we focus on a cavity mode that is resonant with the absorption threshold ('Fermi edge'). 
In agreement with experimental results, the strong light-matter coupling is maintained under very 
high doping within our ladder diagram approximation. While the lower polariton is qualitatively 
unaltered, it acquires a finite lifetime due to relaxation of the valence band hole if the electron 
density exceeds a certain critical value. On the other hand the upper polariton has a finite lifetime 
for all densities, because it lies in the electron-hole continuum where no bound state exists. Our 
calculations show that a narrow upper polariton quasiparticle still exists as a result from the interplay 
between light-matter coupling and final state Coulomb interaction. 



Microcavity polaritons are bosonic quasiparticles that 
appear in a semiconductor microcavity with an embed- 
ded quantum well (QW) when the cavity is tuned to 
an excitonic transition in the QW. In the last decade 
they have emerged as quasiparticles with many favorable 
properties for studying quantum optics and many-body 
physics in integrated photonic structures. Bose-Einstein 
condensation of polaritons and their superfluid properties 
have been investigated in detail both experimentally and 
theoretically [1]. In these studies, the tunability of the 
polariton properties is often exploited. For example, the 
number of quantum wells is varied to change the Rabi fre- 
quency; etching, strain and surface acoustic waves have 
been used to create potentials for the polaritons. 

A tuning parameter that has received relatively little 
attention is the introduction of charges in the QW that 
interact with the polaritons. In the case of low charge 
concentration, trionic bound states exist and a trionic 
polariton is observed 0, An experiment with high 
density modulation doping was performed by Gabbay et 
al. 0| • They have shown that strong light-matter cou- 
pling is possible, even though neither an excitonic nor a 
trionic state could be resolved. 

From the theoretical side, the mixed electronic- 
polaritonic system was proposed to reach supercon- 
ductivity at higher temperatures, possibly even room 
temperature, thanks to a strongly attractive electron- 
electron interaction mediated by the polaritons [5[. For 
what concerns the effect of the electron gas on the polari- 
tons, it has been suggested to facilitate polariton lasing 
thanks to the polariton-electron scattering [6]. In this 
work, the polariton quasiparticle itself was assumed not 
to be affected by the electron gas. A first theoretical 
study on the modification of the single particle polari- 
ton properties due to the electron gas started from an 
estimate of the oscillator strength computed in the con- 
text of the Mahan singularity [I, H[. Indeed, the optical 
excitation of an electron gas (metal or doped semicon- 
ductor) is a central problem in many-body physics that 
was introduced by Mahan Q and to which Nozieres, De 
Dominicis [To| and Anderson [ll[ made seminal contri- 



butions. It is still an active topic of research that has 
recently attracted the attention of the community work- 
ing on ultracold atomic gases 

The purpose of this Letter is to study the different 
regimes of light-matter coupling in highly doped semi- 
conductor QWs, see Fig. [TJ The regime of low doping 
density, where trionic correlations are dominant, is not in 
the scope of this work. A first effect of the presence of the 
high density electron gas is that the low-lying momentum 
states cannot be used for the formation of the exciton, 
because of the Pauli principle. More importantly, the 
electron gas has, in contrast to the empty QW, low ly- 
ing particle-hole excitations. They are responsible for the 
screening of the electron-hole Coulomb interaction. Both 
the Pauli blocking and the screening reduce the exciton 
binding energy 15,[l6j. As long as the particle- hole ex- 
citations are virtual, they do not qualitatively alter the 
light-matter coupling, but merely renormalize the Rabi 
frequency. Qualitative changes occur when the polariton 
energy becomes sufficiently high to excite free electron- 
hole pairs. 

A systematic approach to study the different regimes 
is diagrammatic perturbation theory [17j . In this formal- 
ism, excitonic effects are taken into account by including 
the electron-hole Coulomb interaction through ladder di- 
agrams (see Fig. [2]). The screening of the Coulomb inter- 
action can be taken into account by a dielectric function, 
that we approximate by the RPA result for the jellium 
model. Furthermore, we assume that the quantum well 
electrons and holes are purely two-dimensional. To be 
specific, we will consider n-doping, so that only a free 
electron gas is present and a single photo-created hole. 
The primary physical quantities of interest, the absorp- 
tion and photoluminescence spectrum, can be character- 
ized through the photon spectral function, defined as 
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Figure 1. Microcavity containing a highly n-doped QW. The 
system oscillates coherently between the state with a photon 
in the microcavity and the state with a valence band hole and 
an additional electron in the QW. 



propagator, given by 
D(iuj) 



iuj — uo c — Yl{iuS) 
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Here uj c corresponds to the cavity mode which is assumed 
to be lossless. The photon momentum is taken to be zero 
for simplicity. The photon self energy H(iuj) originates 
from the material excitations within the quantum well 
and thus contains all many-body physics of the electron 
gas. We will approximate it by the Feynman diagrams 
shown in Fig. [2j The Coulomb interaction in these di- 
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agrams is already taken to be the RPA screened inter- 
action, indicated by the double dashed lines. While for 
the electron a complete Dyson series for the screened ex- 
change interaction is taken into account, for the hole both 
the Dyson series for the propagator and the dynamical 
self energy are expanded to lowest order to keep the cal- 
culations manageable: 
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E (k, iv n ) = ^exchange(k) + I!i6u66Je(k, iv n ). (4) 

Obviously, this expansion is only valid for small self en- 
ergy. The approximation (j3j) , (j4]) can be shown to be 
equivalent to the use of Fermi's golden rule to calculate 
the linewidth of a polariton due to the relaxation of the 
hole by the creation of one particle-hole excitation in the 
electron Fermi sea (see supplemental material). The dy- 
namical self energy of the electron is omitted, because its 
relaxation is suppressed by Pauli blocking. 

Under the assumption that the final state Coulomb 
interaction (FSI) is frequency-independent, a Bethe- 
Salpeter equation (BSE) for the electron-hole propagator 
can be derived |9l. It reads 
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Figure 2. (a) Dyson series for the photon propagator D (iuj). 
(b) The photon self energy U(icj) is obtained in the ladder di- 
agram approximation. (c,d) Electron G e and hole Gh Green's 
functions used in the calculations. 
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Here, the Hamiltonian Hx (jSj) describes the relative part 
of the two-particle system; the reduced mass is m r = 

+ m h 1 )- 

To arrive at the BSE, we had to make the static ap- 
proximation for the FSI, for which we used the 2D RPA 
expression of Stern jTsj. This static approximation sim- 
plifies the subsequent calculations, but it will be interest- 
ing to explore the physical consequences of relaxing this 
approximation. A well known limitation of RPA is that 
it is only valid for a high density electron gas a^kp » 1 



with <2q being the effective Bohr radius. Specifically for 
this problem, it means that RPA is not able to capture 
trionic effects that become important at low density. 

The photon self energy is finally related to the vertex 
factor P(k, iuo) as 

5 k 

Here g is the light-matter coupling and S the quantum 
well surface. 
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It can be seen from Eq. (j5j) that the bound excitonic 
levels show up as poles of P(k, iuo) and consequently as 
poles in the photon self energy U(ioj). For an undoped 
Q W this can be calculated analytically [l9[ ; here we im- 
plement it numerically. The lowest lying exciton level is 
depicted in Fig. [3^ as a function of clx^f (blue line), 
where ax /do = (1 + m e /m/ l )/2 holds. Its energy is mea- 
sured in exciton Rydbergs Rx = h 2 /2m r a x . 

A first important energy scale in the exciton physics is 
the lowest energy of a free electron- hole pair, indicated as 
in Fig. [3b and by a red line in panel a. To be consis- 
tent with the computation of the bound exciton, we have 
included the renormalization of the Fermi energy due to 
the exchange interaction. The zero of energy in Fig. [3^i 
was chosen as the energy of a free electron-hole pair with 
both particles at the Fermi wave vector. This energy cor- 
responds to the absorption threshold E t h ('Fermi edge') 
in Fig. [3b. 
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Figure 3. Below axk F the exciton forms the groundstate 
of the optically excited quantum well. Above the critical 
wavevector the exciton becomes metastable (dashed line) and 
obtains a finite lifetime due to relaxation of the hole via emis- 
sion of a particle- hole pair (right figure). axk F is determined 
by rrie/mh. In the figure m e /mh = 0.13 is used. 

Below a critical density, the electron and hole form 
a stable bound state below £? e h, the well-known exci- 
ton. In this case, the hole self energy ([7]) vanishes at the 
bound state energy, corresponding to an infinite exciton 
life time. Beyond the critical value axk Fl which depends 
on the ratio m e /mh (trend indicated with a black arrow), 
the exciton crosses the free electron-hole energy and is no 
longer the lowest energy state, but becomes a metastable 
bound state (blue dashed line). It is the hole kinetic en- 
ergy that is responsible for the fact that the exciton is not 
the ground state. When the hole is bound to the electron 
in a zero momentum exciton, the momentum of the elec- 
tron has to be compensated by an opposite momentum 
of the hole. Because the low momentum electron states 
are Pauli blocked, there is a minimal kinetic energy cost 
k F /2rrih for the hole to form a bound state with the elec- 
tron. Obviously, this cost vanishes in the limit where 
the hole becomes infinitely heavy — »• oo. We then 



find that the exciton is stable for arbitrary density and 
corresponds to the bound state that always exists for an 
attractive potential in 2D. 

When the exciton is metastable, the hole self energy ([7]) 
is nonzero at the exciton energy. It is seen to correspond 
to a relaxation process through the emission of a low en- 
ergetic particle-hole pair, see Fig. [3b. Note that only the 
lifetime contribution of this process is taken into account 
and the associated energy shift has not been considered 
in this work. The delta function in ([JJ expresses energy 
and momentum conservation for this process, while the 
occupation factors describe the amount of phase space 
that is available for the decay process. Below the critical 
Fermi wave vector k F , this phase space factor vanishes 
when energy and momentum are conserved. The relax- 
ation of the hole by plasmon emission was considered 
in the plasmon-pole approximation [20|. However, this 
relaxation channel does not contribute because it does 
not conserve energy and momentum. Unfortunately, the 
hole self energy Q diverges for all energies higher than 
£?th ('the continuum') so that we cannot use our previous 
approach to include the hole self energy. A more sophisti- 
cated approach is necessary. As a first approximation, we 
follow Schmitt-Rink [2l| and obtain a first meaningful re- 
sult by omitting the hole self energy altogether. Actually 
above E t h, the photon self energy acquires an imaginary 
part even when the Hamiltonian (|6]) is real. Physically, 
this finite lifetime comes from the conversion of photons 
into free electrons and holes. Therefore, as a lowest order 
approximation to the spectral function in this regime we 
calculate the spectral function exactly in the same way 
as before, summing up all ladder diagrams to infinite or- 
der, but the self energy for the hole is now discarded. 
The solution of the BSE equation can now be calculated 
analytically to some extent [21[ but we implemented it 
straightforwardly numerically. For the exciton, we re- 
cover the 'excitonic enhancement' [2lL l22|. an increase in 
the absorption close to Eth- This is a many-body mod- 
ification, due to final state interactions, to the constant 
absorption obtained in a single particle picture. 

A summary of the exciton physics is presented in Fig. 
HH- It shows the photon spectral function as a function of 
energy and the electron Fermi wave vector, in the limit 
of very large positive cavity detuning. In this limit, we 
recover the exciton physics. The exciton is blue shifted 
when the electron density increases. It is also clear that 
the excitonic linewidth then increases. The excitonic en- 
hancement for uj > Eth is not very clear on this figure 
and is shown in more detail in the inset on Fig. [4j>e. 
It is an intriguing question whether this excitonic effect 
may lead to the formation of polaritons. Experimentally, 
this question was answered positively in Ref. [4] where 
the observation of 'Fermi edge polaritons' was claimed. 

To answer this question theoretically, we plot in Fig. 
Hb the photon spectral function when the cavity photon 
is in resonance with the Fermi edge E t h- Different detun- 
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Figure 4. a) Photon spectral function as a function of electron Fermi wave vector and energy (arb. units). The free electron-hole 
energy E e h is indicated with the solid white line in the left lower corner, b) shows the photon spectral function at uj c — Eth> The 
lower and upper 'Fermi edge polaritons' are clearly visible. Different cuts along ax^F (white dashed lines) show the narrowing 
of the photon spectral function (black solid) under the continuum (red line) with respect to the exciton spectral function (black 
dashed). The excitonic enhancement is shown in more detail in the insets. Note the different scales showing that the upper 
polariton has a much narrower lineshape as compared to the excitonic enhancement. 



ings are shown in the supplemental material. Several cuts 
at fixed electron densities, indicated with white dashed 
lines, are shown in the panels (c-e). As expected, the 
lower polariton is pushed below the excitonic resonance. 
The shift of the polariton with respect to the exciton is 
however reduced due to the presence of the electron gas. 
This is expected since the light-matter coupling is pro- 
portional to the electron-hole overlap, which is reduced 
due to Pauli blocking and screening. Importantly, the 
linewidth of the polariton is strongly reduced with re- 
spect to the excitonic linewidth. This narrowing can be 
attributed to the reduced final state phase space that 
slows the hole relaxation. For low electron gas densities 
it is even possible for the polaritons to decrease their 
energy below the free electron-hole energy E e \^ making 
their linewidth vanish (on the figure a very small non- 
zero value is added to make the resonance visible). The 
linewidth reduction can be seen more clearly from the 
various cuts taken for different values of clx^f (Fig- IU> 
e). Qualitatively however the lower polariton is entirely 
similar to the one in the empty quantum well. 

This is not the case for the upper polariton that can be 
clearly identified in Fig. HJd, above the absorption thresh- 
old. Above i^th bound states do no longer exist and one 
should think of the excitonic enhancement as the valence 
band hole interacting with all electrons in the Fermi sea 
simultaneously. Nevertheless, this polaritonic resonance 
has a linewidth that is much smaller than the excitonic 
enhancement feature (see insets of Fig. |4j>e). Moreover, 
it is shifted from the cavity resonance by AE ^0.01 Rx • 
Despite the two different approximations that were used 
to compute the spectral function below and above the 
absorption continuum (red vertical line) respectively, the 
matching of the curves is satisfactory. 

It is conceptually interesting to elucidate the role of 



the final state interactions on this upper polariton. To 
this purpose, we have turned off the final state Coulomb 
interaction [Vq = in Eq. flg|)] in our calculations. With- 
out FSI a much smaller energy shift from the cavity reso- 
nance is observed and a larger linewidth is obtained. It is 
really the final state Coulomb interactions between the 
photo-excited electron-hole pair that pushes the quasi- 
particle above the cavity frequency. We have also checked 
that the 'Rabi splitting' vanishes when the cavity energy 
is blue detuned with respect to E t h- It is indeed the 
Fermi edge E t h that is responsible for the formation of 
an upper polariton in the continuum, see supplemental 
material. Finally, it is worth mentioning that when the 
light-matter coupling is increased, the upper polariton 
is further blue shifted and its linewidth increases. For 
a good visibility of the upper 'Fermi edge' polariton, it 
is actually advantageous if the Rabi splitting is not too 
large (not too many quantum wells). This trend is oppo- 
site to the one of the lower polariton, where a larger Rabi 
splitting leads to a lower energy and therefore a longer 
lifetime. 

Note that the upper polariton energy decreases in Fig. 
Ht> for increasing fcp, even though the photon energy is 
constant and the exciton energy increases. In an undoped 
quantum well, the opposite dependence is observed for in- 
creasing exciton energy. We attribute the decrease of the 
upper polariton energy to the reduction of the electron- 
hole overlap when the density is increased. 

Our theoretical results are in qualitative agreement 
with the experiments in Ref. [4], where at high doping 
density both an upper and lower polariton were observed. 
In this experimental work, a phenomenological explana- 
tion based on the Fermi edge enhancement of the optical 
response was suggested. We have shown here that this 
enhancement is due to the usual excitonic bound state for 
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the lower polariton, where it is of a different, collective 
nature for the upper polariton. 

The discontinuity in the spectral function at Eth. is of 
course an artifact of our different approximations below 
and above this continuum energy. It should be cured by 
a consistent inclusion of the self energy below and above 
the absorption threshold energy. This self energy should 
actually contain the diagrams that are relevant for the 
physics of the Mahan singularity so to capture the ef- 
fect of this physics on the nature of the polaritons in 
doped quantum wells. The present work is only a first 
step to such a complete description. It will also be of 
great interest to investigate the interactions between po- 
laritons in doped quantum wells. It may be expected 
that interactions that are induced by excitations in the 
electron gas may enhance the polariton-polariton interac- 
tions with possible applications in quantum information 
processing. 
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EXCITON HAMILTONIAN CONTAINING FINITE LIFETIME FOR THE VALENCE BAND HOLE 

In this section we derive the expression for the matter Hamiltonian (6) used in the article. 

With the assumption of a frequency-independent final state Coulomb interaction, from the ladder diagrams in Fig. 
2, a Bethe-Salpeter equation for the electron-hole bubble P (k, iu) can be found along the lines of Ref. [17]. This 
integral equation for the vertex factor P (k, iuj) is given by 

F- 1 (k, iu) P (k, iu) - ^k- kl P (ki, ^) = g. (9) 

ki 

The solution of this equation yields the photon polarization by the following relation 

n M = f E p ( k > iw )> ( 10 ) 

k 

where for simplicity we have set the photon momentum equal to zero. The function F (k, iu) is given by the Matsubara 
summation over the internal fermionic frequency iv n of the electron and hole propagators, defined in Fig. 2. Because 
we consider a zero photon momentum, only electrons and holes with opposite momentum contribute. Hence we have 

F(k,iu) = -- ^ G e (k, iu - ivn) G h (-k, iv n ) . (11) 
For the electron we use the free Green's function 

Ge(k,«/ n ) = — L— , (12) 

with = + Ea;(k) — fi e . Here, = h 2 k 2 /2m e is the single particle dispersion in the parabolic mass approximation 
and Y< x (k) stands for the screened exchange interaction of the optically excited electron with the Fermi sea. Also the 
electron chemical potential \i e is renormalized due to exchange interaction. 

For the hole Green's function, we include the interactions with the electron Fermi sea perturbatively as follows: 

G h {k,iv n ) = G h (k,iis n ) + G° h (k,iis n )X RPA (k,iis n )G° h (k,w n ). (13) 

Here, for the dynamical self energy Yl RPA , describing hole relaxation that is responsible for a finite polariton lifetime, 
we use the 2D RPA result. The free hole Green's function is given by 

G°(k> n ) = (14) 

The hole dispersion is also assumed to be quadratic, so = h 2 k 2 /2rrih — H>h- No exchange contribution shows up 
because we consider a single hole (fih = 0). 

In order to be able to perform the Matsubara summation in eq. ([TT]) exactly, we further approximate the RPA 
result to lowest order in the Lindhard polarization ('1 particle- hole bubble'), so only emission of a single particle- hole 
pair is taken into account: 

Z RPA (k,iv n ) = Ei6u66Ze(k,2i/ n ) + 0(2 bubbles). (15) 
Using the 2D expression of the Lindhard polarization function, the dynamical self energy is given by 



h 

k-q 



SW(k,^n) = -^E (^) 2 [Mefc) ~ n F (ek^)] "J ^ _T'~" ^ ■ (16) 

° K ,q n ^k-q £ K+q £ K 

Here, the exchange correction to the electron energy was neglected to simplify the subsequent calculations. Further- 
more, the double summation runs over the center-off-mass momentum and relative wavevector of the particle- hole 
pair, resp. q, K. 

Combining the expressions ([TT]) -([T6 ]) we can do the Matsubara sum over iv n . The result is given by 

F(k,icj) = S (l + ^), (17) 
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where we have 



So (k, iu) 



1 - n F (g) - n F (£ k ) 



(18) 



and 



Si (k, iu) 2 ^ f Tr s\ 2 / e w, , e v 

5^— y = ^E(^) ^(^ +q )[l"^( £ K) 



q,K 
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1 1 1 

-^-k"(-^+«) -^k-(-^k+ q +efc+ q -e&) 



RS) [-^ k -(-^ k+q +^ +q --K)] 2 



~ 



i-^fe) -^ k+q +^k +q -^k+^-^k (-^ k+q +^ +q -^+^ k ) 2 J 



(19) 



Upon analytical continuation iu ^ u — E g — [i e — [ih + irj + only the last line yields an imaginary part and thus the 
lifetime of the hole. The other terms correspond to the energy shift, but these are beyond the scope of this work. 
One has 



Im 



Si 2itu- e h _ 



So S* 1 - n F (£) 



n F l-e_k+ q + e k+q- e K 



q,K 



k+q £ K+q + & 



e a h \ 

k £ -k y 



k+q "T- ^K+q — fc K 



(k)] , 



where the delta- function expresses conservation of energy and momentum. This is in fact expression (7) used in 
the article. The 4- fold summation can be worked out further to a 2-dimensional integral which we than compute 
numerically. This leads to the following result (all energies expressed in exciton Rydbergs Rx)' 



lm S ~ Ak) u-k*-X x (k) 



where the Pauli blocking factor S(k) is 



s(k) = e (* > * F ) e (a, > y-^-j^ + s,(k)/jfc 



(20) 



(21) 



Here the hole self energy is 



-'hole 



(k) 
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* d0 ^ F ~ a ~ ~ ^ max (f mi "' fc ^) 2 - ^min Q (a < 0) 9 (fc| - a > iCiJ) (22) 

[9 + »(«)] 



whith s(q) the inverse screening length as defined in The functions a(]<L,q,6),K mm (]!L,q,9) are given by 

a(k, g, 61) = — <? 2 - 2^fcg cos 6» + — (fc 2 - w + S x (k)) 
rrih rrih m r v y 

if min (k,g,(9) = — |a(k,g,0) -<? 2 | 



(23) 



For the BSE the inverse of the function F(k,iu) was needed, see eq. (j9j). Combining eqs. (fT7]h ([18]) and (|20|) we find 

-l 



(u-k 2 -Z x (k)) (l + iS(k) 



= (k) 



w-fc 2 -£ x (k) 



(24) 



For a small hole self energy, we can use the following Taylor expansion: 1/(1 + x) ~ 1 — x for x « 1. Identifying 
x = iE (k) £^ e (k)/ (co> — k 2 — E x (k)) we obtain Fermi's golden rule: 



F-\k,u) = (w-fc 2 -£ x (k)) (l-*S(k) 



J hole 



(k) 



u - k 2 - E x (k) 
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= uj - k 2 - E,(k) - i E (k) Efc oZe (k). (25) 

Plugging this expression into. eq. (j9]) the identification with eq. (6) from the article is straightforward. It is worth 
noting that the Taylor series is consistent with the expansion of the RPA self energy to lowest order in the Lindhard 
polarization. Moreover, the polariton spectrum in figure 4 from the article shows that the linewidths are much smaller 
than the exciton binding energy. Hence, the current results and approximations are selfconsistent. 



INFLUENCE OF FINAL STATE INTERACTIONS AND FERMI EDGE 



Both final state interactions and the presence of the Fermi edge are necessary for the formation of the upper Fermi 
edge polariton. This is illustrated quantitatively in Fig. [5] (a-b) and (c-d) for two electron densities and cavity 
detunings. The blue lines indicate the single-particle results while the red lines contain the final state Coulomb 
interactions. 

Figures a,c show the photon spectral function on resonance with the Fermi edge Et^ (black dashed line) for two 
different electron gas densities. A clear suppression of the spectral function linewidth when taking into account final 
state interactions (FSI) is observed as compared to the single-particle result. 

To estimate the importance of the Fermi edge, the cavity mode is largely blue detuned from the absorption threshold 
energy to E t ^ + 100 (black dashed lines), see figures b, d. It is clear that both the single particle theory and final state 
interaction results now coincide (same scales are used to facilitate comparison). The Fermi edge is thus responsible 
for the 'Rabi splitting', meaning that it really pushes the upper polariton into the electron- hole continuum. 

As a last remark, for very large blue detuning of the cavity mode from the Fermi edge, a red shift from the cavity 
mode is seen in figures b, d. This can be attributed to a different index of refraction in the electron-hole continuum 
in comparison with the excitonic region (see also Ref. [23|). 



a x k F =1.0135 , a x k F =2.5946 
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Figure 5. Comparison of photon spectral functions with (red) and without final state interactions (blue) for different Fermi 
wave vectors. In order to appreciate the importance of the FSI the red and blue curves in figures a-c are to be compared. There 
a significant narrowing of the linewidth can be observed when FSI are taken into account. The role of the Fermi edge for the 
formation of the upper Fermi edge polariton is clear from the comparison of a-b and c-d. The 'Rabi splitting' vanishes for a 
highly blue detuned cavity mode. Both the FSI and the Fermi edge are necessary for the formation of the upper polariton. 



